Fermi arcs in cuprate superconductors: tracking the pseudogap below Tc and above T* 
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Using an energy-momentum dispersion for Bi2Sr2CaCu2 08 obtained from angle-resolved photo- 
electron spectroscopy we show that the shrinking Fermi arc model of the pseudogap is inconsistent 
with Raman scattering below Tc and specific heat near T*. By simulating the quasiparticle energy 
dispersion curves we show that Fermi arcs are an artifact of a T-dependent scattering rate. 

PACS numbers: 74.25.Gz, 74.25.Jb, 74.62.Dh, 74.72.-h 
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The normal-state properties of underdoped high-Tc 
cuprates are dominated by a gap in the density of states 
known as the pseudogap (PG)]l|, Its origin is un- 
known. Unlike a superconducting gap which closes at T^, 
angle-resolved photoemission spectroscopy (ARPES) and 
tunneling in underdoped cuprates have long suggested 
that the PC fills as (1 — T /T*) and disappears abruptly 
at T* [i] . If so then the lost low-energy spectral weight 
would be restored at T* with important thermodynamic 
consequences, as we shall see. The gap energy, Eg, and 
T* decrease with increasing doping, falling to zero at a 
critical doping p = pcrit ~ 0.19 holes per Cu[2]- Gen- 
erally T* is much greater than Tc but when p exceeds 
Pi 0.16 then T* falls below Tc@. 

Recent ARPES studies [1, Q suggest that for T < T* 
the PG covers only part of the Fermi surface (FS) near 
the (7r,0) zone boundary. This leaves ungapped arcs on 
the FS ("Fermi arcs") which grow with increasing T. The 
arcs extend to an angle 6o given by 6'o = f (1 — T/T*), 
where is measured from (7r,0). The gap is thus nodal 
at T = and, with increasing T, retreats towards (7r,0) 
where it closes abruptly at T*. Assuming that the PG 
continues to exist as the underlying normal state below 
Tc, the T-dependent restoration of the pristine FS would 
have important testable consequences, as follows: 

(i) if the PG retreats abruptly to (7r,0) at T* the specific 
heat coefficient, 7, will exhibit an anomaly. The area un- 
der the anomaly is exactly equal to the restored entropy. 
We calculate this below for several Fermi arc scenarios 
and show this is not observed in the experimental data; 

(ii) if the Fermi arcs grow as T rises then the spectral 
weight taken up by the PG falls. So the superfluid den- 
sity would first increase, then fall as T approaches Tc, 
contrary to the observed monotonic decrease [1]; 

(lit) and if the PG closes at T* then in the doping range 
0.16 < p < 0.19, where T* < Tc, the PG should be 
open at T = but closed when T — Tc. Again this 
is not observed, as we show in Fig. [T] Here Tc{p) and 
T*{p) = Eg/2kB are plotted for BiaSraCaCuaOg+a (Bi- 
2212). Panel (a) shows two properties at Tc'. the jump 
in A7 and the electronic entropy S(Tc). Panel (b) shows 
two ground state T = properties: the superfluid density 
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FIG. 1: Thermodynamic data for Bi-2212 (a) two properties 
at Tc: the jump in specific heat coefficient A7 and the elec- 
tronic entropy S{Tc) in mJ/g.at.K; and (b) two properties at 
T = 0: the superfluid density Aq ^ and Xcrit the critical Zn 
concentration giving Tc = 0. Vertical dashed lines indicate 
where T* = Tc and pcrit = 0.19 where the PG closes. 
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the critical density of Zn required to sup- 



press Tc. The data is from Tallon et al. 9]. The two ver- 
tical dashed lines show the region where T* < T^. In (b) 
both ground-state properties abruptly reduce when the 
PG opens at pcrit ■ Now if the PG were to fully close at T* 
then the two T — Tc properties displayed in (a) should 
not begin to fall until T* — Tc i.e. at the vertical line 
at p « 0.16. They do not. In fact, like the ground-state 
properties, they fall abruptly at Pcrit = 0.19. This sug- 
gests a relatively T-independent loss of spectral weight 
due to the PG below Tc, impacting roughly equally on 
ground-state and T — Tc properties. The Fermi arcs, if 
they exist, cannot continue to collapse around the nodes. 

Here we explore these points further. We first con- 
sider the low-T behavior of the Fermi arcs by calculat- 
ing the Raman response. Below Tc the SC gap opens 
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FIG. 2: (Color online) Calculated Raman response at r=0 
for a Fermi arc length fixed below Tc. (a) nodal B2g and (b) 
antinodal B\g symmetry. Arrows show positions of SC and 
PC gap features, (c) PC and Fermi arc given by Eqn. [2] at 
temperatures Ti and T2. (d) eo(T) from Eqns. [3 Hand [5] 



up on the Fermi arcs thereby obscuring the details of 
how the PG further evolves with T, and in particular 
whether the PG itself becomes nodal at T = 0. Raman 
scattering allows the PG and SC gaps to be separately 
probed via the antinodal (-Big) and nodal (-623) Raman 
response 
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1^ . We calculate these for Bi-2212 us- 
ing an ARPES derived energy-momentum dispersion. By 
comparing with Raman data we are able to confirm that 
the PG does not evolve significantly below Tc- 

We then turn to the second issue as to how the PG 
evolves near T* . Does it really close at T* as widely 
believed? We use the same dispersion to calculate the 
specific heat and show that the Fermi arc model leads 
to a large anomaly in 7(T) at T* that is not observed. 
We go on to suggest these issues may be resolved by 
incorporating quasiparticle (QP) lifetime broadening. 

We employ the six parameter tight-binding Bi-2212 
dispersion, e(k), reported by Norman et al.^^ and as- 
sume a rigid single-band approximation. Inclusion of 
band splitting will not significantly affect the fol- 
lowing results. We take the Fermi level to be 34meV 
above the (tt, 0) van Hove singularity (vHs) near optimal 
doping[l3| and 96meV above the vHs in an underdoped 
sample with p = 0.1l[ll|. We interpolate between for 
intermediate doping levels. The imaginary part of the 
unscreened non-resonant Raman response at T = is 
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FIG. 3: Experimental Raman data for Hg-1201 from Le Tacon 
et aZ[l^. (a) -B29 and (b) -Big response for T <^Ta. 



given by[16[ 
Xo (9 = Oi^) 



-i^.(.-2i.(k))^ 



l7(k)|^ 



(1) 

where the integral is over occupied states below Ep, 
A(k) = iAo (cosfca; — cos ky) is the d-wave SC gap func- 
tion and E{k) = ^/e{ky + |A(k)|2. In the Big scatter- 
ing symmetry 7(k)^i9 = 7-Big (cosfc^: — cos ky), giving a 
dominant response from the antinodal sections of the FS. 
For i?2g, 7(k)^^9 = jB2g sin kx sin ky and the response 
is mainly nodal. The magnitude of the SC gap, Aq, 
is taken from the weak-coupling result 2Ao = 4.28kBTc 
and Tc is given by the empirical relationfl7l| Tc/Tc^max = 
1 - 82.6 (p - 0.16)^ We adopt a PG of the form 
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where < 6 < tt/2. We initially assume 
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(3) 



and T* 



E, 



jg^max/ks- This form of the PG is illus- 
trated in Fig. [2Jc) and (d). Eqn. [3] models the lin- 
ear T-degendence of the Fermi arc length inferred from 
ARPES[7]. At T=0, 9o = 7r/4 and the PG is fully nodal. 
As T rises, 6q decreases resulting in a 'fiUing-in' of the 
PG and the growth of the Fermi arcs. The PG is a non- 
states-conserving gapfl^l i.e. unlike the SC gap there is 
no pile up of states outside the gap. This is implemented 
by removing states with e(k) < Eg from the integration 
in Eqn.[T] The doping dependence of Eg is obtained from 
the reported leading-edge ARPES gap at 100X0. 

Figure [2ja) and (b) show the nodal (-829) and antin- 
odal (Big) Raman response for six dopings spanning the 
range 0.12 to 0.19. We consider two scenarios: 
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FIG. 4: Calculated Raman response at T = with a fully 
nodal PG for (a) nodal and (b) antinodal Big symmetry. 
The intensity scale is the same for all plots. 



(i) Firstly, we have assumed that the length of the 
Fermi arc becomes fixed at the onset of superconductiv- 
ity, implemented by setting T = Tc in Eqn. [3l Fig. [2l^a) 
shows the nodal (-B23) response. Leaving aside the 
anomalous electronic Raman continuum above the pair- 
breaking gap, the calculations closely resemble the recent 
results of Le Tacon et o/.[l3] which are reproduced in 
Fig.m The PG peak maximum in the Big response shifts 
monotonically to higher energies with decreasing doping. 
Simultaneously the intensity of this peak rapidly reduces 
with underdoping. In contrast, the SC peak maximum 
in the i?2g response is found to shift to lower energies in 
the underdoped regime. The magnitude of the S23 peak 
persists relatively undiminished down to the lowest dop- 
ing levels. Also reproduced is the increased linear slope 
of the response at very low doping. 

(ii) Secondly, we show in Fig. [4] the Raman response 
in the alternative case where we have assumed that the 
Fermi arcs continue to collapse below Tc. This is done 
setting T = in Eqn. [3] resulting in a fully nodal PG. 
In this case the -623 peak shifts monotonically to higher 
energies with decreasing doping and the intensity reduces 
rapidly. This behavior is not observed experimentally. 
We cannot say that the Fermi arc freezes exactly at Tc 
but our analysis indicates that it remains finite at r=0 
and evolves only weakly below Tc- 

We turn now to the question as to how the PG evolves 
above Tc as T ^ T*. It is widely assumed that the PG 
closes at T* , as indeed Eqn. [3] suggests, thus exposing 
a pristine FS. In fact, Kanigel et a/. [3] suggest that a 
discrete jump in Oq occurs at T* so that the PG retreats 
to the flat sections of the FS near (tt, 0) from whence it 
abruptly disappears at T = T*. This model is given by 



1(1-0.68^) (T<r*) 
(T > T*) 



(4) 







□ Exp. data p=0. 138 

Eqn. 3 

Eqn. 4 

Eqn. 5 



CO 

—5 
E 



1 - 




100 

Temperature (K) 



200 



FIG. 5: Specific heat coefficient calculated for 60 given by 
Eqns. (3), (4) or (5) and compared with the experimental 
data for Bi-2212[ii] where every 20th data point is shown. 



entropy to the bare-band value. This, quite generally, 
will result in a 7 anomaly, the area of which equals the 
restored entropy. The lower is T* the greater is the 
anomaly. We evaluate this here. 

Using the method described previously Q, together 
with the above tight-binding dispersion, we have com- 
puted 7(T) for three cases with p « 0.138. In Fig. [5] we 
compare these with experimental data for Bi-2212. The 
three cases are shown in Fig.[2][d). They are (i) the linear 
behavior described by Eqn. [Sj (ii) the sudden jump in 
inferred by Kanigel et al. and described by Eqn. [U and 
(iii) the smooth evolution of the gap given by[8|: 



1 — tanh 



(5) 



and is illustrated in Fig. ^d). The closure of the PG, 
whether according to Eqn. [3] or Eqn. 21 will restore the 



As can be seen, cases (i) and (ii) lead to substantial 
anomalies in 7 which are clearly not found experimen- 
tally. The experimental data, shown by the data points, 
are from a previous study on the specific heat of Bi- 
2212[3| with the closest corresponding doping state. 
The smooth evolution of the PG given by Eqn. [5] sat- 
isfactorily describes the data and generally discounts the 
possibility that the PG closes abruptly at T* . 

We have so far shown that (i) below Tc the Fermi arcs 
do not continue to shrink notably; and (ii) above Tc the 
Fermi arcs cannot spread out to abruptly form a con- 
nected pristine FS at T* . (iii) Both results are confirmed 
by the data in Fig. [1] in the interesting case where T* lies 
below Tc- These problems seriously prejudice the current 
picture of Fermi arcs but could be resolved as follows by 
invoking a T-dependent scattering rate. 

Firstly, one is led to this view by the observations of 
Norman et al.^^. They modelled the QP peak using 
a QP self energy with a scattering rate or inverse life- 
time Fq which, in underdoped samples, grows with T 
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FIG. 6: Simulated symmetrized ARPES quasiparticle EDCs 
ranging from the antinode to the node = to 45°). The 
broadening is fixed at 0.45i?9(0). The crossover from double 
to single peaks (bold curve) presents an apparent, though 
false, closing of the PG and recovery of Fermi arcs. Inset (a) 
true normalized gap and apparent gap for two broadenings 
0.45x and 0.6xi5g(0)). (b) curve: apparent Fermi arc length 
(FAL) assuming broadening oc T/T*; data: from ref. 0] 



and equals the gap magnitude precisely at T* (see Fig. 
2(b) of ref. This would lead to a smearing out of 

the gap at Ep resulting in a single peak in the disper- 
sion at Ep that looks like a recovered pristine FS. But of 
course the gap is still there, and interestingly is reported 
to be T-independent, just as we concluded above. 

In Fig. [S] we have simulated the symmetrized QP EDCs 
as the sum of two Lorentzians. We use a fixed broadening 
of 0.45£^g(0), corresponding to a particular temperature 
below T* . The /c-dependent gap value, Eg{0) is assumed 
to be d-wave-like with Eg{e) = Eg{0) cos{2e). Fifteen 
angles ranging from 6=0 (the antinode) to 0=45° (the 
node) are shown and the data qualitatively reproduces 
the reported experimental data[3, [l^. These EDCs re- 
veal a crossover from double to single peaks, with the flat 
bold curve lying at the boundary. Within the Fermi arc 
model this would be interpreted as a closure of the PG at 
9 « 30° with a pristine Fermi arc extending from 9 = 30° 
to 45°. But, as shown in inset (a), the true gap docs not 
close until the node at 9 = 45°. The apparent gap, found 
by reading off the peak positions, is also plotted in inset 
(a) and this falls to zero at 9 = 30°. We also show the 
apparent gap for a larger broadening (=0.6£'g(0)) which 
would correspond to a higher temperature, closer to T*. 
Here the apparent gap closes at 23°. 

Inset (b) in Fig. [6] shows the apparent (though ficti- 
tious) Fermi arc length obtained from these simulated 
EDCs when the broadening F = r*(r/T*), where F* is 
the critical value that "closes" the gap at the antinode. 
The arc length exhibits a rapid change at T* just like the 
data of Kanigel et al. '{J] which is also plotted. This rapid 



change arises from the flat part of the d-wave gap when 
6* ^ and does not signal an abrupt recovery of the FS. 

Thus QP lifetime broadening with a d-wave gap ac- 
counts for all apparent Fermi arc features, including the 
abrupt jump in arc length reflected in Eqn.lH We expect 
therefore that the T-variation of the PG is not in Fermi 
arcs but in the scattering rate. Both the PG and the 
SC gap parameters should be replaced by complex terms 
of the form Eg{9) ^ ^^^^ and A(0) -^^P^y 
This naturally leads to a "U-shaped" gap[20j, and in- 
stead of frozen Fermi arcs below Tc the Raman data 
would then insist on a frozen scattering rate. 

In summary, we have used an e(k) dispersion and a 
model for the normal-state PG, both based on ARPES 
results, to show that the shrinking Fermi arc picture is 
inconsistent with Raman data below and thermody- 
namic data near T* . Only by freezing the length of the 
Fermi arcs at or near Tc do we find that the calculations 
mimic the experimental data. This implies that the PG 
does not evolve greatly below T^- We have further shown 
that closure of the PG at T* and the associated recov- 
ery of a pristine FS would lead to a large specific heat 
anomaly that is not observed. We thus question the con- 
cept of Fermi arcs and suggest that they are probably an 
artifact of a T-dependent lifetime scattering rate. 
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